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Abstract 

We give a proof of the recently proposed formula for the dyon spectrum in CHL 
string theories by mapping it to a configuration of Dl and D5-branes and Kaluza-Klein 
monopole. We also give a prescription for computing the degeneracy as a systematic 
expansion in inverse powers of charges. The computation can be formulated as a problem 
of extremizing a duality invariant statistical entropy function whose value at the extremum 
gives the logarithm of the degeneracy. During this analysis we also determine the locations 
of the zeroes and poles of the Siegel modular forms whose inverse give the dyon partition 
function in the CHL models. 
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1 Introduction 

There exists a proposal for the exact degeneracy of dyons in toroidally compactified het- 
erotic string theory [1, 2, 3, 4, 5] and also in toroidally compactified type II string theory [6] 
in four dimensions. These formulae are invariant under the S-duality transformations of 
the theory, and also reproduce the entropy of a dyonic black hole in the limit of large 
charges[2]. Ref. [3] proposed a proof of this formula by first relating this to a five dimen- 
sional D1-D5 system and then counting the degeneracy of this system. 

In [7, 8, 9] this conjecture was generalized to a class of CHL models[10, 11, 12, 13, 14, 
15], obtained by modding out heterotic string theory on T 2 x T 4 by a ffli N transformation 
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that involves 1/N unit of translation along one of the circles of T 2 and a non-trivial action 
on the internal conformal field theory (CFT) describing heterotic string compactification 
on T 4 . The values of iV considered in [7] were N =1,2,3,5,7, with N — 1 representing 
toroidal compactification. The conjectured formulae for the dyon spectrum are invariant 
under the duality symmetries of the theory, and also reproduce the black hole entropy for 
large charges up to first non-leading order. 

The goal of this paper is twofold. First we shall give a 'proof of the conjectured 
formula by relating the dyons in the four dimensional CHL models to a D1-D5 system in 
five dimensional CHL models. Although the basic idea behind our analysis is similar to 
that of [3] , the details are different and even for the N — 1 case our proof does not reduce 
exactly to that given in [3] (see [9] for an attempt to generalize the analysis of [3] to CHL 
models). The other goal will be to develop a systematic procedure for extracting the 
asymptotic behaviour of the degeneracy formula as a series in inverse powers of charges 
and non-perturbative corrections. During this analysis we find that the statistical entropy, 
defined as the logarithm of the degeneracy, is obtained by extremizing a function. We 
call this function the statistical entropy function in analogy with the black hole entropy 
function whose extremization gives the entropy of an extremal black hole. The analogy in 
fact goes further since we show that up to first non-leading order the statistical entropy 
function of the dyons matches the entropy function of extremal black holes carrying the 
same charges. 

Since the analysis involves a lot of technical details, we have organised the paper such 
that only the basic ideas are presented in the text and all technical details are relegated 
to the appendices. In §2 we describe how the dyons in four dimensional CHL models 
can be related to a D1-D5 system in five dimensional CHL models. In §3 we count the 
degeneracy of the D1-D5 system under consideration and show that the result of this 
analysis agreees with the conjecture given in [7, 8]. Finally in §4 we develop a systematic 
procedure for extracting the asymptotic behaviour of the degeneracy in the limit of large 
charges. 

2 Dyons from D1-D5 System 

In this section we shall follow the method of [3] to relate the dyons in a four dimensional 
CHL model to a rotating D1-D5 system in a five dimensional CHL model. The steps 
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leading to this relation are as follows. 

1. Begin with type IIB string theory on K3 x S* 1 with Q 5 D5-branes wrapped on 
K3 x S 1 , Qi Dl-branes wrapped on S 1 , —n units of momentum along S 1 and 
angular momentum Ji and J 2 in two independent planes with Ji + J 2 = ./[16]. 1 
Since a K3-wrapped D5-brane carries —1 unit of Dl-brane charge, this system has 
Qi — Qb units of Dl-brane charge along S 1 . For definiteness we shall choose the 
coordinate along S* 1 such that it has period 2ir. 2 

2. Now place this system at the center of Taub-NUT space with coordinates of the 
Taub-NUT space transverse to K3 x S 1 . The orientation of the Taub-NUT space is 
chosen such that the isometry direction of the Taub-NUT geometry that becomes the 
compact direction S 1 in the asymptotic region coincides with the angular coordinate 
along which the black hole rotates[17]. The new configuration now corresponds 
to a state in type IIB string theory on K3 x S 1 x S 1 with Q 5 D5-branes wrapped 
on K3 x S 1 , (Qi — Q§) units of Dl-brane charge along S 1 , — n units of momentum 
along S 1 , a Kaluza-Klein monopole associated with the compact coordinate S 1 and 
J units of momentum along S 1 . 

3. Make an S-duality transformation on this system to get type IIB string theory on 
K3 xS 1 xS 1 with Q 5 NS5-branes on K3 x S 1 , (Q1 — Q5) units of fundamental string 
winding charge along S 1 , —n units of momentum along S 1 , J units of momentum 
along S 1 , and a Kaluza-Klein monopole associated with S 1 compactification. 

4. Now make an R — > 1/R duality transformation along S* 1 to convert the theory 
into IIA on K3 x S 1 x S 1 with Q5 Kaluza-Klein monopoles associated with S* 1 
compactification, (Qi — Q$) units of fundamental string winding charge along S 1 , 
—n units of momentum along S l , J units of fundamental string winding charge 
along S 1 , and a single NS5-brane wrapped on K3 x S 1 . Here S 1 denotes the dual 
circle of S 1 . 

: If we identify the tangent space group SO (A) of the four transverse spatial dimensions to SU(2)l x 
SU(2)r, then J can be identified with twice the E/(1)l generator of SU(2)l- 

2 Note that in our convention supersymmetry acts on the right-moving sector of the D1-D5 world- 
volume theory, i.e. on modes carrying positive momentum along S 1 . Thus a BPS state will involve only 
excitations involving modes with negative or zero momentum along S 1 . 
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5. Finally using the string-string duality[18, 19, 20, 21, 22] relating heterotic on T 4 to 
type IIA string theory on K3, we can map this system to heterotic string theory on 
T 4 x S 1 x S* 1 with Q 5 Kaluza-Klein monopoles associated with S 1 compactification, 
(Qi — Qb) units of NS5-brane charge along T 4 x S 1 , — n units of momentum along 
S 1 , J units of NS5-brane charge along K3 x S 1 , and a single fundamental string 
wrapped on S 1 . If Q e and Q m denote the electric and magnetic charge vectors 
respectively, then this system has | Q 2 m = (Qi — Q^Q^, \Q\ = n, and Q e -Q m = J- 

Now we can mod out both sides of the duality by a transformation which, in the 
original theory, involves 2n/N translation along S l and some order N transformation g 
acting on K3 that commutes with all the space-time supersymmetry transformations. In 
the final heterotic string theory the transformation acts as 27r/N translation along S 1 
together with some action on the coordinates of T 4 and the 16 left-moving internal 
bosonic coordinates. In order to get a invariant configuration so that we can carry 
out the TLn modding, we need to put periodic boundary conditions on all the branes 
which extend along S 1 , and take N identical copies of all the branes transverse to S* 1 and 
place then at intervals of 2-n/N along S 1 . The latter set includes the five branes along 
K3 x S 1 ] we need to take NJ five branes and place them at intervals of 2ix/N along 
S 1 . After orbifolding the direction along S 1 can be regarded as a circle of radius 1/N, 
and we shall have J five branes per unit period transverse to S 1 . The natural unit of 
momentum along S* 1 is now N, and momentum — n along S* 1 can be regarded as —n/N 
units of momentum. The other charges have the same values as in the parent theory. 
This gives \Q\ = n/N. Q 2 m and Q e ■ Q m remain unchanged from their original values. 
Thus we have 

\Ql = n/N, X -Q 2 m = (Q 1 - Q 5 )Q 5 , Q e -Q m = J- (2.1) 

One point to note is that since we have a single fundamental string along S 1 , a 7L^ 
modding would correspond to requiring that under a shift of 27r/N along S 1 the non- 
compact coordinates of the string transverse to S* 1 come back to their originl values and 
the coordinates of the string along K3 come back to the image of the original location 
under the 'Kn generator g acting on K3. This represents a twisted sector state in the 
orbifold theory. This is consistent with the assertion of [7] that the proposed answer 
for the degeneracy of dyonic states is valid only for states which carry electric charges 
compatible with twisted sector states. 
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This analysis shows that the problem of computing degeneracy of dyons in the four 
dimensional CHL model can be reduced to the problem of computing the degeneracy of a 
rotating D1-D5 system in Taub-NUT space, carrying momentum along the direction com- 
mon to the Dl-brane and the D5-brane. We shall now turn to the problem of computing 
degeneracy of this system. 

3 Counting of States of the Rotating D1-D5 System 
in Taub-NUT Space 

In this section we shall compute the degeneracy of the D1-D5 system described in the 
previous section. For simplicity we shall restrict our analysis to the Q 5 — 1 case, - the 
generalization of this analysis to the Q 5 ^ 1 case has been discussed in appendix D. In 
this case an intuitive picture of the states of the D1-D5 system may be given as follows [24]. 
First of all dynamics of the 1-5 strings produces an effective binding between the Dl and 
D5 which forces the Dl-branes to move in the plane of the D5-brane. Thus the only zero 
modes of the Dl-brane associated with the transverse directions are those along K3. A 
generic state of the system will contain a certain number of isolated Dl-branes in the 
plane of the D5-brane, with the i-th Dl-brane wound Wi times along S 1 and carrying 
momenta — U along S* 1 and ji along S 1 , with w^l^ji £ ^- The BPS condition will also 
require Wi > 0, U > 0. The overall motion of the Dl-D5-system in Taub-NUT space can 
carry some additional momentum — 1 along S* 1 and jo along S 1 . Furthermore the low 
energy dynamics of closed string modes localized near the core of the Taub-NUT space 
can also carry some momentum —l' Q along S 1 . 3 Thus we have the constraint 

J2wi = Qi, l' + l + J2k = n, j + J2ji = J ■ (3-1) 

i i i 

Each of the systems described above must satisfy the boundary condition that under a 

translation along S* 1 by 2-k/N its oscillation modes along K3 get twisted by g. 

Our goal is to compute the degeneracy of the D-brane system described above. In fact 

what we shall really compute is not the degeneracy, but an appropriate index defined as 

follows. First of all the Dl-D5-Kaluza Klein monopole system in type IIA on K3 x S 1 x S 1 

3 Localized low energy excitations around the Taub-NUT space do not carry any momentum along 
S . They can carry fundamental string winding charge along S' 1 [23], but for the configuration which is 
of interest to us, these charges vanish. 
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breaks twelve of the sixteen supercharges of the bulk theory. Thus there will be twelve 
fermionic zero modes associated with the broken supersymmetry generators. Quantization 
of these fermion zero modes produce a 2 6 = 64-fold degeneracy of states. These states form 
a single irreducible multiplet of the supersymmetry algebra known as the intermediate 
multiplet containing equal number of fermionic and bosonic states. This is then tensored 
with the states obtained by quantizing the rest of the degrees of freedom of the theory 
to get the full spectrum. If the state that is tensored with the basic supermultiplet 
is bosonic we shall call this a bosonic supermultiplet and if this state is fermionic we 
shall call this a fermionic supermultiplet. In other words we shall call a supermultiplet 
bosonic (fermionic) if the highest spin state in the supermultiplet is bosonic (fermionic). 
We shall denote by h(Qi,n, J) the number of bosonic supermultiplets minus the number 
of fermionic supermultiplets of the system described above carrying quantum numbers 
(Qi,n,J). By an abuse of notation we shall continue to refer to this number as the 
degeneracy of states. 

We shall first analyze the degeneracy (Ikk{1'q) associated with the low energy dynamics 
of the taub-NUT space, carrying momentum —l' Q along S 1 . Under the duality transfor- 
mation described in section 2 the Taub-NUT space gets mapped to a fundamental string 
wound once along S 1 and carrying momentum —l' Q along S* 1 in heterotic string theory 
on (T 4 x S 1 )/ "E N x S 1 . Since this is a twisted sector state, the problem of computing 
dRKiJ-'o) reduces to the problem of computing the degeneracy of twisted sectors states in 
this string theory. This has been done in appendix A and the result is 

oo 

Y, d KK(l'o)e 27Til ^ = 16e- 2mf> [] {(1 " e 27rinf> y^Ti(l - e 2 ™^)"^} (3.2) 
i' n=l 

The Taub-NUT space breaks eight of the sixteen supersymmetry generators in type IIB 
string theory on K3. This gives rise to eight fermonic zero modes. The factor of 16 in 
(3.2) arises from the quantization of these zero modes. 

Next we turn to the degeneracy g?cm(^o?Jo) associated with the overall motion of 
the Dl-D5-system in Taub-NUT space carrying momentum — 1 along S* 1 and jo along S 1 . 
The low energy theory describing the overall dynamics of the D1-D5 system in Taub-NUT 
space is that of a supersymmetric a-model with Taub-NUT target space, the coordinate 
along S* 1 being identified as the world-sheet a coordinate of this supersymmetric field 
theory. Since S 1 / TL^ has period 2tt/N, the natural unit of momentum along a is N. 
As a result a BPS state of the D-brane system carrying momentum — 1 correponds to 
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Lq = Iq/N, L = O. 4 On the other hand the momentum j along S* 1 is the U(l) charge 
associated with the angular direction of the Taub-NUT space which becomes the compact 
circle asymptotically. The degeneracy of such states has been computed in appendix B, 
and the result is 

E d cM{k, j Q )e***>p+™M = 4e" 2 ™ (1 - e ~ 2 ™)- 2 

'o ,jo 

oo 

J~J |^ e 2mnNp^ ^ ^2mnN p+2mv\^~2 ^ ^2-winN p-2mv^ -2 1 



n=l 



(3.3) 



The D1-D5 system in K3x Taub-NUT space breaks four of the eight supersymmetry gen- 
erators of type IIB string theory on K3 x Taub-NUT. Quantization of the associated 
zero modes gives rise to the factor of 4 in (3.3). The factor of 16 in (3.2) and 4 in (3.3) 
together provides the 64-fold degeneracy of a 1/4 BPS supermultiplet. As pointed out 
earlier, these 64 states contain equal numbers of bosons and fermions. After factoring this 
out we count the rest of the states with weight +1 for bosons and weight —1 for fermions. 

Let us now turn to the computation of the degeneracy n(w, l,j) of a single Dl-brane 
moving inside the D5-brane, carrying winding w, momentum — / and angular momentum 
j. We denote by a the coordinate along the length of the Dl-brane, a being normalized so 
that it coincides with the target space coordinate in which the original S 1 had period 2n. 
Since the Dl-brane carries winding charge w, a changes by 2nw/N as we traverse the whole 
length of the string, regarded as a configuration in the orbifold. Under a — > a + 2irw/N, 
the physical coordinate of the Dl-brane shifts by 2nr/N along S 1 where 

r = w mod iV . (3.4) 

TL N invariance then requires that under a — > a + 2nw/N the location of the Dl-brane 
along K3 gets transformed by g r = g w '. 

We expect the low energy dynamics of this D-brane system to be described by a 
superconformal field theory (SCFT) with target space K3 subject to the above boundary 
condition. Since the Dl-brane is subject to twisted boundary condition with period 

4 Throughout this paper L and Z of a state in a CFT will denote the left and right-moving world- 
sheet momenta with a normalization in which the world-sheet coordinate a has period 2tt. In RR sector 
we implicitly subtract the factors of c/24 (c/24) from L (L ) which arise in going from the sphere to the 
cylinder coordinates. 
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2ttw/N, the natural unit of momentum on the brane is N/w. Thus a total momentum 
—I corresponds to —Iw/N unit of momentum and can be identified with the L — L 
eigenvalue of the state. Since the BPS condition on the Dl-brane corresponds to L = 
we have L = Iw/N. Thus we are looking for a state in the SCFT with 

L = Iw/N, L = . (3.5) 

The bosonic and fermionic excitations on the brane satisfy identical boundary conditions 
in order to preserve space-time supersymmetry. Hence the state belongs to the Ramond- 
Ramond (RR) sector. Furthermore the boundary condition described below eq.(3.4) tells 
us that the state is twisted by g r . Finally, since the total momentum along S 1 is — /, under 
translation by 2ix/N along S 1 this state picks up a phase e ~ 2ml / N , Thus the projection 
operator onto 'ZLn invariant states is given by 

-i N-l 
iV s=0 

Putting these results together we see that the total number of 7L^ invariant bosonic 
minus fermionic states of the single Dl-brane carrying quantum numbers w,l,j is given 
by 

-1 N-l 

l,j) = t;E e^ N Tr RR ~ r (g s (-l) F+F 5 N L ,i w 5j,) . (3.7) 

iV s=0 

where Tr RR ~ r denotes trace over RR sector states twisted by g r , F, F are world-sheet 
fermion numbers (which coincide with the space-time fermion numbers) associated with 
the left and right-moving excitations on the Dl-brane, and J is the U(l) generator asso- 
ciated with an appropriate R-symmetry current of this conformal field theory 5 which can 
be identified with the angular momentum operator[16]. Insertion of (— 1) F in the trace 
automatically projects onto Z = states. 
Let us define [8] 

F^\r, z) ee ^Tr RR ~ r (g s (-l) F+F e 2 - L ^ z ) (3.8) 

Then we have [8] 

F^\r,z) = 1a(t,z), 

5 In this case the full R-symmetry group is 50(4) ~ SU(2)l x SU(2)r and J corresponds to twice 
the U(1)l generator of the SU(2)l group. 
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F^ s \r,z) 
F^ rk) (r,z) 



A(t, z) - — ^— B(t, z) E n {t) for 1 < s < (N - 1) , 



N(N + 1) 



N(N + 1) 



A(t,z) + 



N + l 
2 



£a 1^1 B(r,z), 



N(N + 1) " \ N 

for 1 < r < (N - 1), < k < (N - 1) , 



where 



A(t,z) 



1 2 (t,^) 2 | ^ 3 (r,z) 2 | t? 4 (t,z) ; 



^ 2 (r,0) 2 tf 3 (r,0) 2 tf 4 (r,0) 2 
5(r,^)=r ? (r)^ 1 (r,z) 2 , 



and 



12i 24 



27rmiri2T 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



F( t ' s \t,z) has power series expansion of the form[8] 

F {r ' s \r, z)= ]T c (r ' s) (4n - fy^^z ? 

ng 7L/N,je TL 

for appropriate coefficients c ( - r,s ^(4n — j 2 ). From (3.8), (3.13) it follows that 

^ frwtf (r(-l) F+ ^^ ^^) = c^w/N - j 2 ) . 
Hence (3.7) gives 

JV-l 

n(w, I, ]) = e- 2msl/N c (r ' s) (Alw/N - f) , r = w mod JV . 



(3.13) 



(3.14) 



(3.15) 



s=0 



Using the result for single Dl-brane spectrum, we can now evaluate the degeneracy of 
multiple Dl-branes moving inside the D5-brane. Let dm{W, L, J') denote the degeneracy 
of this system, carrying total Dl-brane charge W = J2i Wi, total momentum — L = —J2ih 
along S* 1 and total angular momentum J' = Then a straightforward combinatoric 

analysis shows that 

£ d m (w,L, j'y^w/N+pL+vj') = Y[ (i_ e ^i(™/N+pi+vj)y n ( w M (316) 



W,L,J' 



u,i,je TL 

w>0,l>0 
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Let us now turn to the full system containing the D5-brane and multiple Dl-branes 
in the Kaluza-Klein monopole background. The combinatoric problem to be solved is the 
following. We have total quantum numbers (Qi,n, J), to be distributed among multiple 
Dl-branes moving inside a D5-brane, the overall D1-D5 system and the Taub-NUT space 
according to the relation (3.1). The degeneracy associated with the dynamics of (multiple) 
Dl-branes inside the D5-brane, carrying quantum numbers (W = J2 w i,L = J2k,J' = 
is given by doiiW, L, J'), the degeneracy associated with the overall dynamics of the 
D1-D5 system carrying quantum numbers Iq, jo is given by dcM(lo,jo) and the degeneracy 
associated with the dynamics of the Taub-NUT space carrying quantum number l' Q is 
given by dcM^'o)- If h(Qi,n,J) denotes the total number of bosonic minus fermionic 
supermultiplets (in the sense described earlier) of the mutiple Dl-brane system carrying 
quantum numbers (Qi,n, J) and if we define 



f(p,a,v)= J2 h(Q 1 ,n,J)e 2ni ^ n+a(Ql - 1)/N+iJ) , (3.17) 

Qi,n,J 



then we get 



f(p,a,v) = —e~ 2 ™' N £ d m (W,L,J')e 2 ^ w / N+ " L+iJ ^ 

l> ' W,L,J< 

E dcM(loJo)e^ ilo ^ ij A \T, d KM)e 2ml 'A . 

\lo,jo J \ 1' J 

(3.18) 

The factor of 1/64 in this equation arises due to the fact that a single 1/4 BPS super- 
multiplet is 64-fold degenerate; thus in order to count the number of supermultiplets we 
need to divide the total number of states by 64. Using (3.2), (3.3), (3.15), (3.16) and the 
relations [8] 

c ( °' 0) (0) = |, c^)(-l) = |, 

c^)(0) = 1 (20 - , c (M(_i) = I, , for S = 1, 2, ... (AT - 1) , 

(3.19) 

one can reduce (3.18) to 

N-l _ y^iV-l -2-itisl/N Jr,s)(Aiy_i2\ 

f(p,&,v) = e - 2wi (P +a / N +'") Y[ TT (1 _ e 2ni(ak'+pl+vj)\ L, a =o W K 3) 

r=0 k'e 7L+jf,i,je 7L 

fc',!>0,j<Ofor k' = l=0 

(3.20) 
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The k' — term in the last expression comes from the terms involving dcM(lo,jo) and 
d-KK^'o)- Comparing the right hand side of this equation with the expression for i>fc given 
in (C.2) we can rewrite (3.18) as 



(i\/N)~ k - 2 

f(p,a,v) = - W"{ . (3.21) 

$fc is a weight k modular form of an appropriate subgroup of the Siegel modular group [8]. 
From (3.17), (3.21) it follows that for n > 0, Q 1 > 1, J > 0, 

h(Qi,n,J) = K f dpdadv— exp [— 2m(pn + a(Qi — l)/N + vJ)\ , 

Jc §h(o. a.v) 



(3.22) 



where 

K = -±(iVN)- k - 2 , (3.23) 
and the integration surface C is the three real dimensional surface: 

Imp = M 2 /N, Ima = M^, Imv = M 3 , 
0<Rep<l, < Re & < N, 0<Rev<l, (3.24) 

Mi, M 2 and M 3 being three large but fixed positive numbers. 

Now recall that from the map described in section 2, we have for this system 

Q 2 e = 2n/N, Q 2 m = 2(Q 1 -Q 5 )Q 5 = 2(Q 1 -l), J = Q e -Q m - (3.25) 

Thus the degeneracy d(Q e , Q m ) of the four dimensional black hole is related to the quantity 
h(Qi,n,J) given via eqs.(3.17), (3.21) as 

= K L dpdadv^—^- exp \-m(pQ 2 e N + &Q 2 JN + 2vQ e ■ Q rn ) . 

JC $ k (p,<T,v) 1 J 

(3.26) 

It has been shown in appendix C that 

$ k (p,a,v)=Q k {<jlN,pN,v). (3.27) 
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Thus denning new variables: 

p = or/N, a = pN, v = v, (3.28) 

we can write 

d(Q e , Qm) = K ( dpdadv ~ } exp \-m(pQ 2 m + dQ\ + 2vQ e ■ Q n 

Jc ®k.(o,<j,v) 1 



(3.29) 



with the integration surface C being the three real dimensional surface: 

Im p = Mi , Im a = M 2 , Im v = M 3 , 
0<Rep<l, 0<Rea<N, 0<Rev<l. (3.30) 

This agrees with the proposal put forward in [7, 8, 9]. 

Although this formula has been derived for Q 5 — 1, a more systematic analysis based 
on mapping the low energy dynamics of the D1-D5 system to a conformal field theory 
with the symmetric product of {Qi—Q^Qz + 1 copies of K3 as the target space[24] can be 
used to show that the final formula holds for general Q$. This has been done in appendix 
D. 



4 Asymptotic Behaviour of the Degeneracy of Dyons 

In this section we shall develop a systematic method for computing the behaviour of 
d(Q ei Q m ) given in (3.29) for large charges. As in [1, 2, 7, 25] we shall try to estimate 
the behaviour of this integral for large charges by deforming the integration contour and 
picking up residues from the poles. This removes one of the three integrals. The remaining 
two integrals will then be performed using a saddle point approximation. At the end of 
this process we need to compare the contributions from different poles and identify the 
one that gives the dominant contribution in the large charge limit. 

The poles of the integrand in (3.29) arise from the zeroes of The locations of the 
zeroes and poles of have been determined in appendix E. According to eqs.(E.16), 
(E.19), has second order zeroes at 

^(vp — v 2 ) + bv + n\G — pmi + m 2 = 0, 

b 2 1 

for mi e N Z£, ni e mi, b e 2 7L + 1, m 2 , n 2 G mini + m 2 n 2 + — = - . 

(4.1) 
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Let us define 

A = n 2 , £ = (m,-mi,^&), y = (p,a,-v), C = m 2 , q = {Q 2 e ,Q 2 m ,Q e • Qm) , 

(4.2) 

and denote by • the SO (2, 1) invariant inner product 

(x\x 2 , x 3 ) ■ (y\ y 2 , y 3 ) = x l y 2 + x 2 y 1 - 2x 3 y 3 . (4.3) 

Then we have 

y 2 = y -y = 2(pa -v 2 ), B ■ y = bv + ma - m x p , (4.4) 

and the first equation of (4.1) may be rewritten as 

±Ay 2 + B-y + C = 0. (4.5) 

Picking up residue at the pole forces us to evaluate the exponent in (3.29) 

-m (pQ 2 m + dQl + 2vQ e ■ Q m ) = -in q ■ y , (4.6) 

at (4.5). To leading approximation the location of the saddle point is now determined by 
extremizing (4.6) with respect to y subject to the condition (4.5). This gives 

q + \(Ay + B)=0, (4.7) 

where A is a lagrange multiplier. (4.5) and (4.7) now give: 



A = W „ 9 g2 o^ » V = -Mt + B). (4.8) 



B 2 -2AC' y A \\ 
Since 

b 2 1 

B 2 -2AC = -2(m l n l + m 2 n 2 + —) = -- (4.9) 
due to the last equation in (4.1), we get 



A = ±V C 27, (4.10) 
Thus at the saddle point the exponential e~ t7rg ' y takes the form: 

E = e -i*q.y = e i*(f?/A\+q.B/A) = e (±-V¥n+^-B)/A ^ 
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Now since q ■ B and A are integers, the second term only gives a phase. Hence 

\E\ = e^v 7 ^ = e *\/QlQ 2 m-(Qe-Qm) 2 /n 2 _ ^ 412 ^ 

Note that we have chosen the sign of the square root so that the sign in the exponent is 
positive since this gives the dominant contribution. 

Eq.(4.12) shows that the leading contribution to the integral comes from the saddle 
point corresponding to n 2 = 1. In this case a p — > p + 1 transformation in (4.1) induces 
ni — * n-i + 1) m 2 — y m 2 — m i, and since rii e 2Z, we can use this symmetry to bring the 
saddle point to ni = 0. On the other hand a a — > a + N transformation in eq.(4.1) induces 
mi — > mi — N, m 2 — > m 2 + niiV. Since mi G iV TL, we can use this transformation to bring 
mi to 0. Finally the v — > tJ+l transformation in (4.1) induces b — > 6 — 2, m 2 — > m 2 + 6— 1. 
Since 6 e 2 22 + 1, we can use this transformation to set 6 = 1. m 2 is now determined to 
be zero from the last equation in (4.1). Thus we have 

m i = m 2 — rii — 0, n 2 = 1, 6 = 1. (4-13) 

The corresponding zero of 5>fc is at 

op -v 2 +v = 0. (4.14) 

The exponent of eq.(4.12) for n 2 = 1 gives the leading term in the expression for the 
statistical entropy, but in order to find a systematic expansion of the entropy in inverse 
power of charges, we need to carefully evaluate the complete contribution from the pole 
at (4.14). For this it will be useful to define new variable p, a, v through the relations: 

pa -v 2 pa - (v - I) 2 pa -v 2 +v . . 

P = -^ , o=- -= -, v = ~ = , 4.15 

a a a 

or equivalently, 

v 2 - pa „ 1 _ v-p 

P=o ' a= o ' v= o • ( 416 ) 

2v — p — a 2v — p — a 2v — p — a 

In these variables (4.14) gets mapped to v — 0, and we have, near v = 0[7], 

Mp,a,v) = A7r 2 (2v-p-a) k v 2 f^(p)f^(a) + 0(v 4 ), 

/W(r) = r)(T) k+2 r)(NT) k+2 . (4.17) 
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Also it follows from (4.15), (4.16) that 

dpdadv = {2v — p — c)~ 3 dpdadv . 
Thus the contribution to (3.29) from the pole at pa — v 2 + v = is given by 

d(Q e ,Q m ) ~ ^ dpdadv v- 2 (2v-p-a)- k - 3 (f< k \p)j< k \a)Y 1 



(4.18) 



47T 2 Jo 

exp 



~ Z7r 1 O ^rn + o Ve + o Ve " Qr, 

[2v — p — a 2,v — p — a 2v — p — a 



(4.19) 



where the integration contour C now encloses the divisor v = 0. The correction to this 
formula involves contribution from other poles for which n 2 7^ 1, and are suppressed by 
powers of e~ Q2 . Evaluating the v integral in (4.19) by Cauchy's formula, we get 



-2(Jfe + 3) + 2vri 



dpda 
{p~Taf 
pa 



p + a 



Ql-—_Ql + p - a 



p + a 



exp 



1 



27T 



-^—Q 2 - 



p + a 

Qe ' Qrn 



Qe ' Qr 



p + a p + a p + a 

In f (k \p) - In f (k \a) - (k + 2) ln(p + a) 



Let us now introduce new complex variables a and S through the relations: 

p = a + iS, a = —a + iS . 



(4.20) 



(4.21) 



Then (4.20) may be rewritten as 

d(Q e ,Q m ) * / 

47T J 

7T 



dSda 



2(k + 3) + - {(a 2 + S 2 ^ + Q 2 - 2ag e • Q m ) 

exp ^{(a 2 + 5 2 )g 2 t + g 2 -2og e -g m } 

- In f ik) {a + iS) - In f (k) {-a + iS) - (Jfe + 2) ln(2iS) 



(4.22) 
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So far we have not specified the contour for a and S integration, except that it must 
pass through the saddle point. In the leading approximation the saddle point, obtained 
by extremizing 

^ {(a 2 + S 2 )Ql + Ql- 2aQ e ■ Q m } (4.23) 

occurs for real values of a and S. Thus we can take the contour of integration to be along 
the real a and S axis. If we now define 



t = -a + iS = Ti + ir 2 , 
then (4.22) may be reexpressed as 



(4.24) 



d(Q e ,Q m ) ~ Koi-lf'^ f ^ 



7T 



2(A; + 3) + -|Q e + rg m | 2 

T~2 



cxp 



7T 



2r 9 



|0e + rQ m \ 2 - ln/ (fc) (r) - ln/( fe )(-f) - (k + 2) ln(2r 2 ) 



where 



■1)' 



1 



(4.25) 



(4.26) 



4vr v " / 47riV( fc + 4 )/ 2 ' 

Note that the degeneracy factor is positive or negative depending on whether Q e ■ Q m is 

positive or negative. This is natural from the point of view of a black hole solution since 

a classical black hole is expected to be bosonic (fermionic) for Q e ■ Q m even (odd). 

Identifying \d(Q e ,Q m )\ with e S3tat ^ e '^ m ^ where S sta t denotes the statistical entropy, 

we can rewrite (4.25) as 



S s tat {Qe >Qm) 



T 2 



-F(r) 



(4.27) 



where f — (ri, r 2 ) or (r, r) depending on the basis we choose to use, and 



F(r) 



IT 

2^ 



\Q e + rQ m \ 2 - ln/W(r) - In f^(-f) - (k + 2) ln(2r 2 ) 



+ In j^o (2(k + 3) + ^-\Q e + rQ m \ 2 ^j 



(4.28) 



Note that F(r) also depends on the charge vectors Q e , Q m , but we have not explicitly 
displayed these in its argument. In (4.27) the integration measure d 2 r/ (r 2 ) 2 as well as the 
integrand e F ^ are manifestly invariant under the Ti(N) transformation: 

ar + b 



Qe OjQe bQmi Qr, 



-cQ e + dQ m , T 



CT + d ' 

a,b,c,de 7L, ad - be = 1, a,d—l mod N, c = mod iV . (4.29) 
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Thus S s tat{Qe,Qm) computed from (4.27) is invariant under Ti(N). 

We shall now describe a systematic procedure for evaluating S sta t as an expansion in 
inverse powers of the charges. For this we introduce the generating function: 

e W(T) = I fl e -F(r) + J.^ (4 3Q) 

J r 2 

for a two dmensional vector J, and define T(u) as the Legendre transform of W(J): 



dW{J) 

U) — J ■ U — VV \ J ) , Ui — 

It follows from (4.31) that 



T(u) = J.u-W(J), Ut = ^a. (4.31) 



dT(u 

■J% = 

As a result if 



J( = ^. (4,2) 



^M. = at u = u , (4.33) 

OUi 

then it follows from (4.30)-(4.32), (4.27) that 

r(«o) = -W(J = 0) = -S stat . (4.34) 

Thus the computation of S sta t can be done by first calculating T(u) and then evaluating it 
at its extremum. T(u) in turn can be calculated by regarding this as a sum of one particle 
irreducible (1PI) Feynman diagrams in the zero dimensional field theory with action 
F(f) + 21nT2. Since S sta t is given by the value of the function —T(u) at its extremum, 
we can identify —T(u) as the entropy function for the statistical entropy in analogy with 
the corresponding result for black hole entropy[26, 27]. 

A convenient method of calculating T(u) is the so called background field method. For 
this we choose some arbitrary base point tb and define 

,W B (r B ,J) _ f -F(r B +ff)+.f-ff 



e W B (T B ,J) _ / - ■/ c -b\T B +f]) + J-fi (4.35) 
•/ ( 'm + V2) 2 

r B (r B ,x) = J-X-W B (T B ,J), X i = dWB ^f ,J) - (4-36) 
By shifting the integration variable in (4.35) to r = t b + ff it follows easily that 

W b (t b , J) = W(J) -T B -J, (4.37) 
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and hence 



T B (r B ,x) =r(f B + x). 



(4.38) 



Thus the computation of T(u) reduces to the computation of Tb{u,x — 0)- The latter 
in turn can be computed as the sum of 1PI vacuum diagrams in the 0-dimensional field 
theory with action F{u + ff) + 2 ln(-U2 + V2), with ff regarded as fundamental fields, and u 
regarded as some fixed background. 

While this gives a definition of the statistical entropy function whose extremization 
leads to the statistical entropy, the entropy function constructed this way is not manifestly 
duality invariant. This is due to the fact that since the duality transformation has a non- 
linear action on (ti,t 2 ), the generating function W(J) defined in (4.30) and hence also 
the effective action T(u) defined in (4.31) does not have manifest duality symmetries. Of 
course the statistical entropy obtained by extremizing T(u) will be duality invariant since 
this is given in terms of the manifestly duality invariant integral (4.28). In appendix F 
we have described a slightly different construction based on Riemann normal coordinates 
which yields a manifestly duality invariant statistical entropy function. The result of this 
analysis is that instead of using the function — Tb{t) as the statistical entropy function we 
can use a different manifestly duality invariant function —Tb{t) as the statistical entropy 
function. Tb{tb) is defined as the sum of 1PI vacuum diagrams computed from the action 



where D T , D f are duality invariant covariant derivatives defined recursively through the 
relation: 



for any arbitrary ordering of D T and D f in D™D™F(t). During this computation the 
components (£, £) or (£1,^2) of £ are to be regarded as the zero dimensional quantum 
fields and f B is to be taken as a fixed base point. The result of this computation expresses 
^b(jb) in terms of manifestly duality invariant quantity F{f) and its duality invariant 
covariant derivatives. 

It has also been shown in appendix F that explicit evaluation of f b(t) gives 



-T B (r) = \Q e + T Q m \ 2 -\nf k \r)-\nf ik \-f)-{k + 2) ln(2r 2 ) + constant + 0(Q- 2 ) . 




D t (D™D?F(t)) 
D f (D™D n f F(r)) 



{d T -im/T 2 )(D™D n f F(T)), 
(d f + tn/T 2 )(D™D n f F(r)), 



(4.40) 



2r 2 



(4.41) 
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Up to an additive constant this agrees with the black hole entropy function for CHL models 
calculated in [27, 28, 29] using the quantum effective action of these theories[30, 31]. 
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A Counting Twisted Sector Elementary Heterotic String 
States in CHL Models 

In this appendix we shall compute 

g(r) = J2dKK(l' )e 2 ^ (A.l) 

where cIkk(Io) is the degeneracy of the BPS states of a Kaluza-Klein monopole associated 
with the compact direction S 1 carrying momentum — 1' along S 1 . As discussed in section 
3, under a duality transformation (Ikk{1'q) gets mapped to the degeneracy of BPS states 
of a fundamental string wound once along S* 1 and carrying momentum —1' along S* 1 in 
heterotic string theory on (T 4 xS 1 )/ TL^ x S 1 . Since we have a singly wound string along 
S 1 , after orbifolding it becomes a twisted sector state. For computing g{r) we can now 
proceed as in [32]. In heterotic string theory on T 6 the internal momenta belong to the 
Narain lattice embedded in a 28-dimensional Lorentzian space of signature (6,22) [33, 34]. 
Let us denote by gu the generator of the TL^ transformation and by V\\ the subspace 
of the 28-dimensional momentum space which transform non-trivially under g^. The 
dimension of V\\ is equal to the number of U(l) gauge fields which are projected out by 
the orbifolding procedure. This in turn is given by [7] 

24 

20-2A;, k = — 2. (A.2) 

N+l 

Let L'q denote the contribution to the L eigenvalue of a state from all the left-handed 
bosonic oscillators and also from the components of the momentum along Vm, and L' 
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denote the contribution to the L Q eigenvalue from the right-handed bosonic and fermionic 
oscillators. 6 Then we have 

L -L = L' -L' -l' /N, (A.3) 

since —l' /N is the contribution to L — L from components of the internal momentum 
invariant under g^. Now level matching condition tells us that L — L = 0, whereas BPS 
condition requires that we have L' = 0. Thus on these states 

L'o = l'o/N . (A.4) 

This allows us to reexpress g{r) as 

1 N-l 

9(r) = E Tr' gH (g s H e^ NL o) , (A.5) 

where Tr' gH denotes trace over states twisted by gu carrying all possible momenta along V\\ 
and involving arbitrary excitation of left-moving bosonic oscillators. We can simplify this 
expression by noting that in the twisted sector of the orbifold theory the TL^ projection 
implements level matching. Since we have already implemented the level matching con- 
dition on the states via eq.(A.4), gu acts trivially on these states. In particular the part 
of gu that corresponds to 2ix/N translation along S 1 gives a phase of e - 2ml 'o/ N = e~ 2mL o, 
and this cancels the phase coming from the part of gu that acts on the twisted bosonic 
oscillators. Thus we can rewrite (A.5) as 

g(r/N) = Tr' gH (e 2 ^«) . (A.6) 

Since the L' Q in the exponent comes from oscillator contribution and components of mo- 
mentum along V\\, only the part of gu twisting which involves the 16 left-moving coordi- 
nates and the coordinates of T 4 are relevant for this computation. 
Some useful properties of g(r/N) are listed below: 

1. From (A.6) we see that g(r/N) can be identified as an appropriate partition function 

with spin structure (0, 1) where the spin structure (s,r) represents twisting by g r H 

6 A11 the components of the right-handed momentum are invariant under gn and hence they do not 
contribute to L' . 



21 



and an insertion of g s H in the trace. Since under a modular transformation r 
(ar + b)/ (cr + d) the spin structure (s, r) gets transformed to 



s \ fa b\ { s 
r ) \c d ) \ r 



(A.7) 



we see that the subgroup of SL(2, %) that preserves the set of spin structures ^ 
mod N is determined by the requirement 

a, 6, c, d G 7L, ad-bc=l, b = mod N, d = 1 mod iV . (A. 8) 

This describes the group T 1 (A^). Thus we expect that g{r/N) will transform as a 
modular form of an appropriate weight under the group r 1 (A^). 

2. The weight of the modular form may be determined as follows. First of all note 
that g(r/N) receives contribution from the 24 left-moving bosonic oscillators. This 
contributes —12 to the weight of the modular form. On the other hand g(r/N) also 
involves a sum over a (20 — 2k) dimensional momentum lattice embedded in V\\. 
This gives a contribution of (10 — k) to the modular weight. Thus the net modular 
weight of g(r/N) is 

24 

10 — Jfe — 12 = — fc — 2 = — — . (A.9) 

N +1 v ' 

3. Since g{r/N) receives contribution from only bosonic oscillators, each term in the 
expansion of g(r/N) as a power series in e 2mT / N has positive coefficients and hence 
g(r/N) does not vanish in the interior of the upper half plane. 

4. For large Im(r), g{r/N) behaves as 

g(r/N) ~ \Qe~ 2mT ' N , (A. 10) 

where the exponent —27iir/N comes from the L' eigenvalue of —1/N associated 
with the ground state of the (twisted) bosonic oscillators [7], and the factor of 16 
counts the 8 Ramond and 8 Neveu-Schwarz sector states associated with the broken 
supersymmetry generators. 

5. For small r the behaviour of g(r/N) is controlled by the large Im{r) behaviour 
of the partition function with spin structure (1,0), i.e. by the trace of gHe 2mL '°/ T 
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over the untwisted sector states. Since the untwisted sector ground state has L' Q 
eigenvalue —1, this gives: 

g{r/N) ~ e 2m/T for r -> . (A. 11) 

From this it follows that gfY/TV)^ 1 is a modular form of T 1 (A^) of weight (k + 2) 
without any pole in the interior of the upper half plane, and behaves as 

(gir/N))- 1 ~ -^e 2mr/N , for large Im r, 
16 

~ e- 2wi/r for small r. (A. 12) 

Since g(r/N)~ l vanishes at the cusps r = n + ioo, n e 2Z, and r = 0, we conclude that it 
is a cusp form of T 1 (A^) of weight (k + 2). This leads to a unique choice for g(r/N)[35]: 

1 6 

9(r/N) = -j^j^ = l^T/Ny^r)-^ . (A.13) 

This gives 

oo 

g( T ) = 16e- 2wiT ]l{(l-e 2winT )-^(l-e 2mnNT y^Ti} . (A.14) 

n=l 

For N — 1 this reproduces the standard result g(r) = 16?](r) _24 for the partition function 
of half BPS states in toroidally compactified heterotic string theory [36]. For N = 2 this 
result agrees with that of [37] obtained by explicit computation. 

B Counting States in the Supersymmetric Field The- 
ory with Taub-NUT Target Space 

In this appendix we will analyze the overall motion of the D1-D5 system in the Taub-NUT 
(Kaluza-Klein monopole) space and count the number of states of this system which carry 
momenta — Iq along S* 1 and jo along S 1 . The Dl-D5-system wrapped on K3 x S* 1 in flat 
transverse space has four bosonic zero modes labelling the transverse coordinates and eight 
fermionic zero modes associated with the breaking of eight out of sixteen supersymmetries 
of type IIA string theory on K3 x S 1 . Thus when the transverse space is Taub-NUT, we 
expect the low energy dynamics of this system to be described by a (1+1) dimensional 
supersymmetric field theory with four bosonic and eight fermionic coordinates, with the 
bosonic coordinates taking value in the Taub-NUT target space. Since the world-sheet 
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coordinate a of this field theory is identified with the coordinate along S 1 , and since 
S 1 / TLfq has periodicity 2n/N, the natural unit of momentum along a is N . As a result a 
state of the D-brane system carrying momentum — 1 corresponds to a state in this field 
theory with L — L = Iq/N. 

We begin our analysis of this field theory by writing down the metric of the Taub-NUT 
space: 

ds 2 = ^l + ^(dr 2 + r 2 {d9 2 + sm 2 9d<p 2 ))+R 2 (l + ^ 1 (2# + cos^0) 2 (B.l) 

with the identifications: 

(9, 0, V) = (2tt - 9, + tt, ^ + |) = (9, + 2vr, 0> + vr) = (0, 0, ^ + 2vr) . (B.2) 

Close to the origin the metric reduces to that of flat space R 4 written in terms of Euler 
angles 9,<f),ip and r, while for large r it is that of R 3 x S 1 , with S 1 parametrized by the 
angular coordinate -0. For later use it will be useful to introduce the cartesian coordinates: 

1 r- ( , <A 2 r- 9 ( , <A 

x = 20" cos - cos I -0 + — I , x = 2yr cos - sm I -0 + — I , 

x 3 = 2\fr sin - cos ^0 — , x 4 = 2\fr sin - sin ^0 — ^ . (B.3) 

The metric (B.l) has a global U(1)l x SU(2)r symmetry. The SU(2)r symmetry 
refers to the usual rotation group of three dimensional space and will not be of interest 
to us. The U(1)l symmetry acts as 

+ (B.4) 

with no action on any of the other coordinates. From the point of view of an asymptotic 
observer this is just a translation along the compact circle S 1 parametrized by -0, and 
the corresponding conserved charge is the quantum number jo- On the other hand using 
(B.3) we see that near the origin the ?fj translation acts as simultaneous rotation along 
the 1-2 and 3-4 planes. Thus near the origin the contribution to the j charge can be 
identified as the sum of the angular momentum in the 1-2 and 3-4 planes. 

Let us now analyze the transformation laws of the fermion fields under U(1)l transfor- 
mation. Since the fermions transform under the tangent space group, we need to find the 
compensating tangent space rotation that accompanies the global U(1)l transformation. 
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Let S0(4) T = SU(2)1 x SU(2)% be the tangent space rotation group of the Taub-NUT 
space, - this needs to be distinguished from the global symmetry group described earlier. 
The Taub-NUT space is known to have SU(2) holonomy. We shall choose our convention 
such that the holonomy is in the SU(2)\. To see what this means we use the fact that 
for a suitable choice of the tangent space basis vectors (e°, e 1 , e 2 , e 3 ) the spin connection 
uj ab = u ab dx^ associated with the metric (B.l) is given by[38] 

R\ , „n — I no { 



u 01 = 



+ -) (r + R)- 1 a x , u; 02 = -(r + -) (r + R) 'ay, 



uj 03 = ~\{r + R)~ 2 a z , c 23 = -f (r + R)' 1 o x , 

c 31 = -|(r + R)- 1 a y , c 12 = (^(r + Ry 2 - l] a, , (B.5) 



where a x , a y and a z are a set of one- forms 

a x = cos (2tp) d9 + sin (2-0) sin , 
o y = — sin (2-0) <i# + cos (2*0) sin 9d<p , 

cr 2 = cos9d(j) + 2d^. (B.6) 
Thus near the origin r = we have 

, ,01 _ , ,23 _ 1 , ,02 _ , ,31 _ 1 _ , ,03 _ , ,12 _ 1 m 7 \ 

u; = a; = — <r x , to = to = — <j y , to = to = — <j z . [d. t ) 
2 2 2 

If S afc denotes the generator of the tangent space group then we see from (B.7) that near 
the origin the holonomy group is generated by the combinations S 01 + X 23 , S 02 + S 31 
and S 03 + S 12 . These generate an SU(2) group, - by our convention this describes the 
tangent space S77(2);[ group near the origin. The precise relation between the S a6 's and 
the generators T 1 , T 2 , T 3 of SU{2) T L are 

T i = I( S oi + S 2 3)) T 2 = i(S 02 + S 31 ), T 3 = i(S 03 + S 12 ), (B.8) 
2 2 2 

so that we have 

u ab Y, ab = -o x T x - GyT 2 - a z T 3 . (B.9) 
Using (B.6) we can express this as 

oj ab T, ab = -iUodUn- 1 , (B.10) 
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where 

Now let us consider the effect of a global U(1)l rotation — > + e. Using (B.10), 
(B.ll) we see that this induces a transformation 

Since T 3 is a generator of S'C/ (2) J we see that in order to preserve the action a global U (1) l 
rotation must be accompanied by a compensating SU (2) J transformation. In particular 
the fermion fields x must transform as 

X^e^x, (B.13) 

where T 3 is to be taken in the representation of SU (2) J in which the fermions transform. 

Now we know that the fermions transform in the (1, 2) + (2, 1) representation of the 
tangent space SU(2)J j x SU(2)^ i group. Thus half of the fermions are neutral under 
SU(2)J j and hence also under the global U(1)l- On the other hand since the holonomy is 
in SU(2) r £, there is no coupling of these fermions to the spin connection and they behave 
as free fermions. The other half of the fermions, carrying non-trivial SU(2)J j quantum 
numbers, are interacting and do transform under the global U(1)l group. Furthermore, 
since type IIB string theory on K3 is a chiral theory, the world-sheet chirality of these 
fermions is correlated with the chirality under the tangent space rotation group. In 
particular, the SU(2)J j singlet free fermions are left-moving on the Dl-D5-brane world- 
sheet, and the SU(2)J i singlet interacting fermions are right-moving on the world-sheet. 7 

This shows that the resulting (1+1) dimensional world-volume theory on the Dl- 
D5 system is described by a set of four free left-moving U(1)l invariant fermion fields, 
together with an interacting theory of four bosons and four right-moving U(1)l non- 
invariant fermions. Let us first calculate the contribution to the partition function due 
to the free left-moving fermions. Since these fermons do not carry any jo charge, their 
contribution is given by: 

y if\ _ rp // i\F+F 2mpl +2-irivj \ 

^freeVAV = J- 1 free left —moving fermions VV / /' 

7 This can be seen as follows. Since the bosons are interacting, the free fermions have no bosonic 
superpartner, and hence they do not transform under the unbroken supersymmetry of the D1-D5 system in 
K3xTaub-NUT space. Thus they must be left-moving. On the other hand the SU(2)] i singlet interacting 
fermions are the superpartners of the interacting bosonic fields and hence must be right-moving on the 
world-sheet. 
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= 4l](l-e 2lM i 4 , (B.14) 

n=l 

taking into account the identification / = (L — L )/N. The factor of 4 comes from the 
quantization of the free fermion zero modes. The latter in turn can be interpreted as due 
to the four broken supersymmetries of the Dl-D5-system on fC3x Taub-NUT space. 

Now we turn to the interacting part of the theory. Since we are computing an index we 
can assume that it does not depend on continuously varying parameters. Let us take the 
size R of the Taub-NUT space to be large so that the metric is almost flat and in a local 
region of the Taub-NUT space the world-volume theory of the D1-D5 system is almost 
free. In this case we should be able to compute the contribution to the non-zero mode 
bosonic and fermonic oscillators by placing the D1-D5 system at any point in Taub-NUT 
space, - say at the origin. The contribution from the zero modes however is sensitive to 
the global geometry of the Taub-NUT space and should be computed separately. 

Since supersymmetry acts on the right-moving bosons and fermions, in order to get 
a BPS state the right-moving bosonic and fermionic oscillators must be in their ground 
state. Thus as far as the contribution due to the non-zero mode oscillators are concerned, 
we only need to examine the effect of left-moving bosonic oscillators carrying L = lo/N, 
L = and angular momentum j . From (B.3) we see that a translation of ip induces 
simultaneous rotation in the x 1 — x 2 and x 3 — x A plane. Hence we need to compute the 
partition function of free left-moving bosons carrying total angular momentum j in the 
x 1 — x 2 and x 3 — x 4 plane. This is easily done using the result of [39] . The answer is 

Z osc (p,v) = Tr oscillators ((-l)^e 2 ^ +2 ^) 

OO 

~~ II ^ _ e 2mnNp+2mvy^l _ e 2mnNp-2mv^2 ' (B.15) 

The result can be understood as follows. The bosonic oscillators in the 1 — 2 plane and the 
3 — 4 plane can be split into complex conjuate pairs which carry opposite units of angular 
momentum. Therefore two of the bosonic oscillators carry +1 unit of angular momentum 
and the other two carry —1 unit of angular momentum. The partition function weighted 
with the angular momentum is then given by (B.15). 

Finally we turn to the contribution Z z to the partition function from the bosonic and 
fermionic zero modes of the interacting part of the theory. Since the intercting theory has 
four bosonic and four fermionic fields, the dynamics is that of a superparticle with four 
bosonic and four fermionic coordinates moving in Taub-NUT space. Under the holonomy 
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group SU (2)£ of the Taub-NUT space both the bosons and the fermions transform in a 
pair of spinor representations. This system is described by an N = 4 supersymmetric 
quantum mechanics. Thus in order to look for BPS states of the D1-D5 system we need 
to look for supersymmetric ground states of this quantum mechanics. 

So far we have been working at a special point in the moduli space of the CHL string 
theory where the circles S l and S 1 are orthonormal in the asymptotic geometry. In this 
case the BPS mass of the Dl-D5-Kaluza-Klein monopole system is equal to the sum of 
the BPS masses of the D1-D5 system and the Kaluza-Klein monopole. As a result there 
is no potential term in the D1-D5 world-volume action and analysis of bound states is 
difficult. But this is not a generic situation. Once we switch on a component of the metric 
that mixes S 1 and S 1 the BPS mass of the Dl-D5-Kaluza-Klein monopole system is less 
than the sum of the BPS masses of the D1-D5 system and the Kaluza-Klein monopole. 
In this case the D1-D5 system in the presence of the monopole has a potential and the 
system is easier to analyze. On the other hand the analysis of the dynamics of non-zero 
modes will not be affected by this modification since we are computing an index, - hence 
our results for Zf ree and Z osc should remain unchanged. 

The mixing between S* 1 and S 1 can be achieved by replacing the dip term in the 
expression for the metric given in (B.l) by dip + ady where y is the coordinate along S 1 
and a is a small deformation parameter. This clearly remains a solution of the equations 
of motion but gives an r dependent contribution to the yy component of the metric: 



As a result the tension of the D1-D5 system, being proportional to y/g yy , acquires an 
r-dependent contribution proportional to 



to first order in a 2 . Supersymmetrization of this term gives rise to other fermionic terms. 

Thus we have to analyze the dynamics of a superparticle with M = 4 supersymmetry 
moving in Taub-NUT space under a potential proportional to (B.l 7). This is precisely 
the problem analyzed in [40, 41] in a different context. The result of this analysis can be 
summarized as follows. Depending on the sign of the deformation parameter a we have 
supersymmetric bound states for jo > or jo < 0, where jo is the momentum conjugate 
to the coordinate ip. In the weak coupling limit the number of bound states for a given 




(B.16) 




(B.17) 
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value of jo is given by |jo|. If for definiteness we choose the sign of a such that we get 
bound states for positive jo, then this gives the zero mode partition function 

oo 2-wiv 

Z z (v) = Tr zcromodcs ((-lf +V^>) = £ jo = n _' 2 • (B.18) 

30=1 \ l e I 

Since this is invariant under v — > — v we shall get the same answer if we had chosen to 
work with the opposite sign of a that produces bound states with negative j . 

Finally putting all the ingredients together the partition function of states associated 
with the centre of mass motion of the D1-D5 system in Taub-NUT space is given by 

E dc M (lo,Jo)e 2ml ^ +2m ^ = Z ilcc (p)Z sc(p,v)Z z (v) (B.19) 

lo ,30 

= 4e _27 ™ i (l — g- 27 ™)- 2 



oo 

X 



e 2irinNp^4^ ^2mnN p+2mv^ -2 ^ ^2irinN p-2mv\^-2^ 



n=l 



C Proof of 3>jfe(T, U, V) = $ k (NU, T/N, V) 

The Siegel modular form which appears in the expression for the dyon spectrum of the 
Win CHL model is given by [8] 



T, V) = -{i^N)- k - 2 exp [2m (±T + U + V 



N-l ( \\YH=o e- Mls ' N T c^)(Uk>-b 2 ) 

II II ll-exp(2m(k'T + IU + bV))\ 

r=0 I, be 7L,k'e 7L+jj I J 
k',l,b>0 

N-l ( \ \Y. N s=o e27raa ' N c{r ' S) ^ lk '- b2 ) 



II II { 1 exp(2rri(k'T + IU + I A ' ) ) 



r=0 i,6e a.fc'e S-ft 

fc',i,f>>0 



(CI) 



where (k r , I, b) > means fc' > 0, / > 0, 6 G S or fc' = 0, 1 > 0, 6 G Z or fc' = 0, / = 0, b < 
0. The coefficients c < - r ' s - ) have been defined in (3.9)-(3.13). From this definition it follows 
that c^ r ' s ^(4n — b 2 ) = c < -~ r ~ s \4n — b 2 ). Hence the two products in (C.l) give identical 
results and we may rewrite this as 

$ k (U,T,V) = -{i^fN)- k - 2 exp(2m(^T + U + V 



29 



II II \ 1 - exp(27ri(Jfe'T + IU + 6V)) I 



r=0 !,&€ K.fc'e 7/. ■ 

fe',i,6>0 



(C.2) 



For N = 1 this reduces to the familiar weight 10 cusp form of the modular group of genus 
two Riemann surfaces [42, 43, 44]. 

If we define j = Nk' in (C.2), then in the above product r = j mod N, and we get 

MU,T,V) = -{iVN)- k ~ 2 exp(27n(^T + C/ + y)) 

1 - exp(2m{jT/N + IU + bV)) \ 

) 



j,l,b>0 



(C.3) 



Substituting U — > T/N, T — > UN in this product and exchanging the dummy indices j 
and I we get: 

$ k (NU, T/N, V) = -(i^fN)- k ~ 2 exp (o,m (^ T + U + V )) 
IJ 1 1 - exp(2m{jT/N + IU + bV)) \ 



j,l,b>0 



(C.4) 



Now using (3.9) one can check explicitly that 8 

N-l N-l 



e^ ils/N F^ s \t, z)=J2 e- 2m3S ' N F^ s \t, z) , (C.5) 

s=0 s=0 

and hence from (3.13) 

N-l N-l 

e -2mis/N c (j,«)( 4n _ 6 2) = e ~ 2m ^ N c^ s \An - b 2 ) . (C.6) 

Eq.(C.6) shows that the right hand sides of (C.3) and (C.4) are identical. Hence 

MT,U,V) =$ k (NU,T/N,V). (C.7) 



8 For this we need to use the relation J2k=o En {^r) = NEn(t). This follows from the Fourier 
expansion of En(t) given in (3.12). 
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D The D1-D5 System and SCFT with Target Space 

S P K3/ 7L N 



In this appendix we evaluate the degeneracy of the low lying states of the D1-D5 system 
described in sections 2 and 3 for Q 5 > 1. Before the 7L N orbifolding the dynamics of 
the relative motion of Q$ D5-branes wrapped on K3 x S* 1 and Qi Dl-branes wrapped on 
S 1 is captured by the M = (4, 4) superconformal cr-model with the symmetric product of 
Q${Qi — Qb) + 1 copies of K3 as the target space [24]. We shall denote this target space 
by S P K3 = (K3) p /Sp, where Sp refers to the permutation group of P elements. The 
world-sheet coordinate a of this conformal field theory is identified with the coordinate 
along S 1 . We shall first review various aspects of the superconformal field theory with 
target space (K3) p /Sp, and then discuss the effect of the projection that is required 
to describe a Dl-D5-brane configuration on the CHL orbifold. 

Let g be an element of Sp and [g] denote the conjugacy class of g. Then the Hilbert 
space of the SCFT with target space (K3) p / S P decomposes into a direct sum of twisted 
sectors labelled by the conjugacy classes of S P : 



where (w) denotes cyclic permutation of w elements and P w is the number of copies of 
(w) in g. Thus these conjugacy classes are characterized by partitions P w of P such that 



If we denote by H w the Hilbert space of states twisted by the generator u of the 7L W 
group of cyclic permutation of w elements, and projected by the same TL W group, then 
(D.4) shows that for the conjugacy class [g] given in (D.2) 




(D.3) 



w 



The centralizer C g of the conjugacy class [g] given in (D.2) is given by 

C g = S Pl x (Sp 2 x ft? ) x • • • x (S Pa x TL P /) . 



(D.4) 



nf^ = ® w>0 s p -h 



■w 



(D.5) 
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Consider first the Hilbert space H w . This twisted sector is represented by the Hilbert 
space of the sigma model of w coordinate fields Xi(a) G K3 with the cyclic boundary 
condition 

X i (a + 27i)=coX i (a)=X t+1 (a), ie(l,...,w), (D.6) 

where u acts by u : Xi — > Xi + \. Therefore the w coordinate fields can be glued together 
as a single field but in the interval < o < 2nw, moving in the target space K3. Thus we 
now have a string of length 2irw, - commonly known as the long string, - moving in K3. 
The natural unit of momentum on this string is 1/w, and hence if a state of this string 
carries physical momentum —I along S 1 , it corresponds to Lq — Lq eigenvalue Iw. 

At this stage it is worth pointing out that whereas for Q5 = 1 the quantum number w 
can be identified with the winding charge of the D-string, this is not so for Q$ > 1. Thus 
we should not regard the long string as a D-string, - rather it provides some effective 
description of the dynamics. 

Once we know the spectrum of 7i w , the full spectrum of the CFT of the D1-D5 system 
is obtained by taking the direct product of the spectrum of Ti^s and then carrying out 
appropriate symmetrization described in (D.5). 

We now turn to the effect of the 'ZLjy projection that is required in order to get a state 
of the D1-D5 system in the TL^ CHL model. For this we need to pick a 'ZLn invariant 
state from each of the 7i w . Thus we need to first study the effect of the ULn projection 
in each 7i w and then take the direct product of these 7L^ invariant subspaces followed 
by appropriate symmetrization. In particular if we denote by 7i' w the subspace of 
invariant states in H w , then the full 7L N invariant Hilbert space of the D1-D5 system will 
be given by 

® W>G S P -H' W . (D.7) 

The E N projection on a single long string is implemented as follows. Since the long 
string has length 2nw, we need to impose periodic boundary condition on this string with 
period 2nw/N. Since under a translation by 2nw/N on the string the physical coordinate 
along S* 1 gets shifted by 2nr/N where r = w mod N, under such a shift the coordinates 
of the string along K3 must be transformed by g r where g is the generator of the 7L N 
action on K3. Thus if ^at denotes the group generated by g, we can view the long string 
as a closed string of length 2nw/N with target space K3/ E N , belonging to a sector of 
this SCFT twisted by g r . Since the natural unit of momentum along the string is now 
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N/w, a physical momentum — / along S 1 would correspond to L — L eigenvalue of Iw/N 
in the SCFT describing the dynamics of the long string. 

Now from (3.14) we know that for an SCFT with target space K3/ 7Ln 

\ Tr RR ~r(g s (-l) F+P 5 NLo , w 5j,) = c^(Alw/N - f) . (D.8) 

Here J is the angular momentum operator. Also the projection operator for TL n invariant 
states with physical momentum —I along S* 1 is: 

s 

Hence the total number of bosonic minus fermionic states in the single long string Hilbert 
space, carrying momentum — / along S 1 and angular momentum j is given by: 

1 J- e-^ ls Tr RR ~rg s (-l) F+ ^ NUh aj,) = E e-^ ls c^\4lw/N - f) ee n(w, l,j) . 

s s 

(D.10) 

According to (D.7) the next step is the evaluation of the partition function for the 
symmetrized tensor products of the Hilbert spaces 7i' w . For this we use the following 
formula from [45]. If d sym (P w ,w, L, J) denotes the number of bosonic minus fermionic 
states in S Pw Ti,' w carrying total momentum — L along S* 1 and total angular momentum J, 
then 

E E^ (P W L J^ e 27riL P+ 27riJi >+2TiyP w _ j~J ^ _ e 2my+2mlp+2-irijv^ n W>JJ 
P w =0 L,J TL 

l>0 

(D.ll) 

The proof of this given in [45] holds provided that the coeficients n(w,l,j) are integers. 
Therefore to apply this formula for the Hilbert space H' w we need 

N-l 

n(w,l,j) = ]T c^ s \Alw/N -j 2 )e-^ ils ' N (D.12) 

s=0 

to be an integer. This can be explicity verified using the formulae given in (3.9)-(3.12). 

Using the identity in (D.ll) we can evaluate the generating function for the bosonic 
minus fermionic states for the relative dynamics of the D1-D5 system. Eq.(D.7) shows 
that all we need to do is to take the product over w of the right hand side of (D.ll). 
More specifically, if h'(P, n, J) denotes the total number of bosonic minus fermionic states 
carrying total string length 2nP/N = 2n J2 w >o wP w /N (counting a single string with 



33 



quantum number w to have length 2nw/N), total momentum — n along S 1 and total 
angular momentum J, then we have 



J2 h\P,n,J)e 2m ^ n ^ p - l ^ N+i,J ^ 

J~J T~J ^ _ e 2iriaw/N+2TTilp+2iTijvy 



P,n,J 



»>() i>0 
JV-1 

= e -27ria/iV -Q -Q ^ _ e 27ri(p'a+ip+jt))^-X; s=0 e^^'V'^) (4p'i-j 2 ) _ 

r=0 p'e m+ft.l.je 7L 

p'>0,l>0 

(D.13) 

We have multiplied e - 2ma l N on both sides of the above equation to absorb the shift of 
unity in the value of P. Also in arriving at the last expression in (D.13) we have defined 
p' = w/N. Physically the quantum number P corresponds to P = Qs(Qi — Q5) + 1. 

Let us now evaluate the full degeneracy h(P, n, J) of the D1-D5 system and the Taub- 
NUT space, counting each supermultiplet as one state. To do this we need to put in the 
contribution from the excitation modes of the Kaluza-Klein monopole and the center of 
mass motion of the D1-D5 system. These have already been computed in appendices A, 
B. Putting together all the results we get 



f(p,a,v) = h (P,n,J)e 2 ^ n+ ^ p -^ N+dJ 



P,n,J 

1 N ~ 1 w_i 

J_ e -2nia/N -Q Yl _ e 2^i(p'a+lp+jv)yY^ N s J e- 2 ^ l / N c^(A P 'l~f) 

p'>O,l>0 

E dcM(lo,Jo)e 2m(lo ^A \ Y,dKM)e 2ml 'A (D.14) 

\h,jo J \l'o J 

The factor of 1/64 in this expression accounts for the fact that each supermultiplet con- 
tains 64 states. Using (3.2), (3.3) one can reduce the above equation to 

N-l 

ffa ^ = e -2ni(a/N+p+v) TT TT ^ _ e 2™{p' B+lp+jvy^J e- 2 ™ l / N c^ (4p'l-j 2 ) 

r=0 P '= r Z+j T ,ie 7L,j 

p'>0,i>0 

(D.15) 

where the product over j runs over integer values for (p', I) 7^ (0, 0) and runs over only 
positive integers (or negative integers) when p' — I — 0. Note that the p' = terms in the 
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above equation arise from the terms involving cIcm{Io, jo) and cIkr^)- Now comparing 
the right hand side of this equation with (C.2) we can write 



f(P>*>*) = ~ ¥p^ (D.16) 

From this point onward the analyis of the degeneracy of the D1-D5 system in Taub-NUT 
space proceeds as in section 3 with Q 2 m = 2Q 5 (Qi — Q 5 ). We have thus extended the 
analysis of section 3 to the Q 5 > 1 case. 

E Zeroes and Poles of <E>£ 

As has been shown in [8] the Siegel modular form <5> k (p,a,v) given in (C.2) satisfies the 
relation: 



-2 In $ fc - 2 In $ fe - 2 k In det ImVt + constant 
1 r d 2 r 

r,s=0l=0 ^ 2 miim2i n 2 e a,ni6 S+J,5e2 K+! 



b 2 

2niT(m 1 n 1 + m 2 n 2 + —) 

4 



x 



where 



and 



exp \n 2 (ap-v 2 ) + bv + n x d - pm x + m 2 f^ e 2nimis/N h\ r ' s \r) , 

(E.l) 

fi=(~ ~) , F = det/mfi, (E.2) 

The functions F^ r ' s \r, z) have been defined in eqs.(3.9). The coefficients c^ r,s \An) intro- 
duced in (3.13) are related to h\ r ' s \r) through the expansion: 

4 r ' s) (r)= c {r ' s \4n)q n , h ( {' 8 \r) = ]T c (r > s \4n)q n , q = e 27riT . (E.4) 
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From eqs.(3.9)-(3.13) it follows that the only non-zero c^ r ' s ^(4n) with negative n are 
c (M(_i) f or a n n, and c (r ' s) (-l + for N — 5,7. They have the values: 

c (M (-l) = 

c (^)(_l) = for r ^ mod N , 

c(M(_i + l) = o, 

c^(-l + |) = - N{N + 7 ){N -l) e2mk/N > forr^OmodiV. (E.5) 

Eq.(E.l) shows that the zeroes and poles of appear only when the r integral on 
the right hand side of this equation diverges from the region near r = ioo. Now, if we 
consider a term proportional to e 2nmT in the expansion of h!f' s \ then for large r 2 , the T\ 
integral gives a non- vanishing answer only if 

n + m\n\ + m 2 n 2 + — = . (E-6) 

Thus after performing the T\ integral, the only r 2 dependence of the integrand in the large 
r 2 region comes from the 

7TT 2 

factor in the exponent. Since this is negative, the only way the integral can diverge from 
the large r 2 region is if this vanishes: 

n 2 (<rp — v 2 ) + bv + n{a — pm\ + m 2 = , (E.8) 

for some mi, m 2 , ni, n 2 , 6 appearing in the sum in (E.l). 
Now we have the identity 

b 2 1 

mini + m 2 n 2 + — = -(p 2 L - p 2 R ) , (E.9) 



n 2 (ap — f ) + 6f + riicr — pm x + m 2 (E.7) 



where 



niiop — v 2 ) + bv + ni<r — pmi + m 2 



— |m 2 + n 2 (a 2 p 2 + <7ipi - t>i - t> 2 ) - mipi + ni0i + k>i} 



If /~~ ~~ ~~ 2^ \ _ / 2v 2 \ V 

+ 2Y I n2 1 ~~ aiP2 + 2 ~ — j + mif>2 + Hl \ ° 2 ~ ~ I ~ I 

+2^ + m?-^i + n 2 ^) 2 . (E.10) 
[2 p 2 p 2 J 
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Since p 2 L is positive semi-definite, and since p\ vanishes when (E.8) holds, (E.9) shows 
that we must have 

m\ri\ + m 2 n 2 + — > . (E-H) 

Furthermore the equality sign holds only when p\ also vanishes. This requires m\ = m 2 = 
ni = n 2 = b = 0. The corresponding divergence is present for all a, p, v and is removed 
by a subtraction term[8]. Thus the divergences which depend on a, p, v come from those 
values of rrii, rii, b which satisfy (E.8) and for which 

mini + m 2 n 2 + — > . (E.12) 
4 

This, together with eq.(E.6), now show that we must have 

n<0. (E.13) 

In other words the only terms in the expansion of hf'^ responsible for a divergent contri- 
bution to the integral (E.l) are the ones involving negative powers of q. From eqs.(E.4), 
(E.5) it now follows that for all N there is a divergent contribution to (E.l) from the 
n = —1/4 term, i.e. the term 

c (r ' s) (-i)<T 1/4 , (E.14) 

in hi' s \ On the other hand for N = 5, 7 there are additional divergences from n = — \ + 
i.e. the term 

c^(-l + ^)g" 1/4+1/Ar , (E.15) 

in h^K Since the subscript of h is odd in both cases, we must have b odd. 
First consider the contribution from the c ( - r ' s \— 1) term. (E.6) now gives 

b 2 1 . . 

mini + m 2 n 2 + — = - ■ (E.16) 
4 4 

After estimating the r 2 integral in (E.l) for n 2 (ap — v 2 ) + bv + nio — pmi + m 2 ~ 0, one 
easily finds that the divergent contribution is given by 

N-l 

-2 e 2nimis/N c (r ' s) (-l) In n 2 (ap-v 2 )+bv + ma-pmi + m 2 



s=0 



r = niN mod N, b = 1 mod 2 , (E.17) 
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where we have included a factor of 2 due to the fact that the lattice vectors (m, n, b) and 
(— m, —ft, —6) give identical divergent contribution. Thus near this region ^ behaves as 

„ / xV^ -1 e ^ im \ s l N c ( r - s )(-l) 

~ (n 2 (ap-v 2 ) + bv + n 1 a-pm 1 + m 2 ) s=0 . (E.18) 

Using (E.5) we can compute the sum in the exponent. In particular r must vanish for 
c^ r ' s \— 1) to be non-zero, and hence n\ must be an integer. Furthermore, substituting 
c(°' s )(— 1) = 2/N in (E.18) we see that the sum in the exponent vanishes unless mi is an 
integer multiple of N. The final result after performing the sum is 

$ fe ~ (n 2 (ap — v 2 ) + bv + n\a — pm\ + m 2 ) 

for mi G N 2£, n\ G K, 6 G 2 S + 1, 777.2, «2 £ ^ , 
~ 1 otherwise. (E.19) 

For N = 5, 7 we also have divergent contribution to (E.l) from the c^ r ' s \— 1 + -^) term. 
In this case (E.6) gives 

6 2 1 1 

777i?7i + m 2 n 2 + — = - - — . (E.20) 
The divergent contribution takes the form 



N-l A 

-2 e 2nimis/N c {r ' s) (-l + — ) In n 2 (ap-v 2 ) + bv + n 1 a-pm 1 + m 2 

s=0 ^ 



r = n x N mod N, b = 1 mod 2. (E.21) 

Thus behaves as 

$k ~ (n 2 (tp - v 2 ) + bd + ni a - p mi + m 2 )^ 3=0 N . (E.22) 
We can compute the sum in the exponent by rewriting it as 

N-l A N-l A 

J2 e 2*im lS /N c (r,s)^_ 1 + _) = ^ ^im.rk/N c (r,rfc) + x for r ^ mod iV . (E.23) 

s=0 ^ fc=0 ^ 



Using (E.5) we now get 

e 2wim lS /N c (r,s)^_i + _ J -48/(iV 2 - 1) for = -1 mod iV (E 24) 

N 1 otherwise 
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This gives 

/ o „ \-48/(N 2 -l) 

$k ~ [n 2 {crp - v ) + bv + n X (j - pm l + m 2 J 

for mimN = -1 mod N, mN ^ mod N, b G 2 H + 1, mi, m 2 , n 2 G S , 
~ 1 otherwise. (E.25) 

To summarize, for iV = 1,2, 3, 5, 7, has a second order zero at (E.8) for mi = 
mod N, rii = mod 1 and odd b, satisfying (E.16). On the other hand for N = 5 and 
N = 7,<& k also has poles of order 48/(iV 2 - 1) at (E.8) for m^N = -1 mod N, niN ^ 
mod N and odd b, satisfying (E.20). (Note that both for N = 5 and N = 7, 48/{N 2 - 1) 
is an integer, and hence the singularities are poles and not branch points.) 

F Riemann Normal Coordinates and Duality Invari- 
ant Statistical Entropy Function 

In section 4 we considered ff = f — tb for some fixed base point tb as the fundamental 
field in defining Wb(t b , J) and T b (tb, x)- 111 this appendix we shall try to generalize this 
by treating 

i=M (F-i) 

as a fundamental field where g(ff) is an arbitrary function of ff with a Taylor series expan- 
sion starting with the linear terms (i.e. g(0) = 0). In this case the generating function of 

— # 

£ correlation functions will be given by 

e W B (T B ,J) = I d2 V e - F (T B +ff)+J-g(rj) _ /p_ 2 ) 

J (TB2 + V2) 2 

Thus Wb(tb,0) = S s tat- The corresponding effective action is 

r B (?B,j) = J-j-WB(TB,I), A = dW f^ J) ■ (F.3) 

OJi 



From this definition it follows that 



Now suppose is a specific value of tb for which 



J, = . (F.4) 



ip=0 



= 0. (F.5) 

j=o 
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In this case we have J = for tfj = 0, and hence 

f B (f(°\ 0) = -W B ^°\ 0) = -S stat . (F.6) 

We shall now show that Tb(tb,0), regarded as a function of tb, has an extremum at 
t b = fW>. From (F.3), (F.4) we see that 

f B (^ 0) + e, 0) = -W B (t* 0) + e,J = df B (r^ + e, $)/df) r ff . (F.7) 



Now 



e ^ s (r s +?,J) = f d2 V e - F (T B +e+rf)+J-g(rj) = f d 2 1] ^p^+^+j.^^ 

J {jB2 + 62 + V2) 2 J (TB2 + m) 2 

(F.8) 

where in the second step we have made a change of variables ff — > 17 — e. Since gifj—e) = 
g(fj) + 0(e), this shows that 

Wb(t b + e, J) = W B (r B , J) + 0(e,J fc ) . (F.9) 
Using this information in (F.7) we get 

(F.10) 

Eq.(F.5) shows that the second term on the right hand side of this equation is of order 
e 2 , and J appearing in the argument of the first term is of order e. Expanding the first 
term in a Taylor series expansion in J, and noting that the linear term vanishes due to 
(F.5), we get 

f b (t*V + e, 0) = -W b {t^\ J = 0) + 0(e 2 ) = f B (^°\ 0) + 0(e 2 ) . (F.ll) 

Thus 

dT f^ =Q at r = r^. (F.12) 

on 

Using (F.6) and (F.12) we see that the statistical entropy is given by the value of —Tb{t, 0) 
at its extremum f = f^°\ Thus we can identify — Tb(t, 0) as the statistical entropy 
function. This is computed as the sum of 1PI vacuum amplitudes in the theory with £j 
regarded as the fundamental fields. 
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We shall now show that for a suitable choice of the coordinates £, the statistical entropy 
function — Tb(t, 0) defined this way can be made manifestly duality invariant. This is 
done by choosing £ as Riemann normal coordinates. For a given base point tb we define 

— * 

the coordinate £ for a given point r in the upper half plane as follows. We introduce the 
duality invariant metric on the upper half plane 

ds 2 = {r 2 )- 2 {drl + dr 2 ) , (F.13) 

and draw a geodesic connecting tb and f. The coordinate £ corresponding to the point 
t is then defined by identifying |£| as the distance between tb and f along the geodesic 
and the direction of £ is taken to be the direction of the tangent vector to the geodesic at 
Tb- 9 Since the metric (F.13) is invariant under a duality transformation, it is clear that if 
a duality transformation maps fg to f' B and f to r', then the Riemann normal coordinate 
£' of f' with respect to t' b will have the property that = |£|. Thus £ and £' are related 
by a rotation. In order to determine the angle of rotation, we note that under a duality 
transformation (4.29), 

dr' = (ct + d)- 2 dr. (F.14) 

Thus 

dr' \ct + d\ 2 dr 



\dr'\ (cr + d) 2 \dr\ ' ^ F ' 15 ^ 
This shows that a geodesic passing through tb gets rotated by a phase \ctb + cI\ 2 / (ctb + cI) 2 
under a duality transformation. Hence 

e-r^^^t' < F - 16 > 

(ct b + dy cr B + a 

where 

£ = 6 + *6, ? = £ + (F-17) 

Since for given tb the duality transformation acts linearly on £, the corresponding 
generating function Wb(t b , J) and the effective action f b (tb, tp) will be manifestly duality 

— * — # 

invariant under simultaneous transformation of tb, J or ip and of course the charges Q e 
and <5m- In particular the 1PI vacuum amplitude Tb(t, 0) will be duality invariant under 
the transformation (4.29). 



9 Often one uses the convention that the distance along the geodesic is y/ gij{fB)€, l €,^ ■ This definition 
differs from the one used here by a multiplicative factor of T2B- Since this transforms covariantly under 
a duality transformation, both choices of £ would give manifestly covariant Feynman rules. 
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We shall now give an algorithm for explicitly generating duality covariant vertices in 
this O-dimensional field theory. For this we need to expand the duality invariant function 

— # 

F(f) in a Taylor series expansion in £. This is given by: 

OO 1 

F (r) = E -rfa*) B fci • • • & D h ■ ■ ■ D tn F(r) , (F.18) 



n=0 n[ 



T = T B 



where Di denotes covariant derivative with respect to Tj, computed using the afline con- 
nection T l jk constructed from the metric (F.13). We arrive at (F.18) by using the result 

— # 

that in the £ coordinate system symmetrized covariant derivatives are equal to ordinary 
derivatives. Using this we can replace ordinary derivatives in the Taylor series expansion 
to covariant derivatives with respect to £j. In the second step we use the tensorial trans- 
formation properties of covariant derivatives to convert covariant derivative with respect 
to & to covariant derivative with respect to Tj. The (r B2 ) n factor in (F.18) arises due to 
the fact that near r = tb, 

dr i = T B2 di i . (F.19) 
In the (r, f) coordinate system the non-zero components of the connection are 

r; r = f , T*„ = -f • (F.20) 
The curvature tensor computed from this connection has the form 

R l jk i = -(fan - 5\g 3k ) , (F.21) 

which shows that the metric (F.13) describes a constant negative curvature metric. From 
(F.20) it follows that 

D T {D™D n f F{r)) = (d T -im/T 2 )(D?D?F(r)), 

D f {D™D n f F{r)) = (d f + in/T 2 )(D?D?F(r)) : (F.22) 

for any arbitrary ordering of D T and D f in D" 1 D™F(t). (F.22) provides us with explicit 
expressions for the covariant derivatives of F appearing in (F.18). Also using (F.22) one 
can prove iteratively that under a duality transformation 

(cr + d\ m ' n 

(r 2 ) m+n D™D n f F(r) -> ——- (r 2 ) m+n D™D n f F(r) , (F.23) 

\CT + dl 
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under a duality transformation. This establishes manifest duality covariance of the ver- 
tices constructed from (F.18). 

We also need to worry about the contribution from the integration measure. The 
original measure rf 2 r/(r 2 ) 2 was duality invariant. Since duality transformation acts on 
as a rotation, d 2 ^ is also a duality invariant measure. Thus we must have 



d 2 T 

(r 2 ) 2 



= J(r B ,Od 2 ^ 



(F.24) 



for some duality invariant function J{t B -,£,)- ^ nas Deen shown in appendix G that 

1 



J(r B ,0 = ^smh\^\ 



(F.25) 



— # 

We can now regard — \uJ{t b ,^) as an additional contribution to the action and expand 

— * 

this in a power series expansion in £ to generate additional vertices. Using the expression 
for F(f) given in (4.28) we now see that the full 'action' is given by 



F(r)-\nJ(r B ,0 



\Qe + rQ m \ 2 - ln/W(r) - In f^(-f) - (k + 2) ln(2r 2 ) 



7l 



+ \n\K -\Q e + rQ r 



7~2 



+ In 1 + 



2(A; + 3)r 2 
n\Q e + rQ m \ 2 



(F.26) 



In this expression the first term inside the square bracket is quadratic in the charges, the 
last term contains terms of order Q~ 2n for n > 1, and the other terms are of order Q°. 
Thus in order to carry out a systematic expansion in powers of inverse charges we need to 
regard the first term as the tree level contribution, the last term as two and higher loop 
contributions and the other terms as the 1-loop contribution. We can now evaluate the 
effective action Tb{tb) in a systematic loop expansion. The leading term in the effective 
action is then just the first term in (F.26) evaluated at f = f B : 



r (r B ) 



7T 



2t. 



\Q e + r B Q r 



(F.27) 



2B 



At the next order we need to include the tree level contribution from the rest of the terms 
in the action (except the last term which is higher order) and one loop contribution from 
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the first term. The former corresponds to these terms being evaluated at f — tb, i.e. 

— # — # 

£ = 0. Since J{f B ,i = 0) = 1, we get this contribution to be 

In /W(r B ) + ln/ W (-f B ) + (fc + 2) ln(2r 2B ) - In \k —\Q e + r B Q m \ 2 )) . (F.28) 

For the one loop contribution from the first term in the action we need to expand this 
term to quadratic order in £ using eqs.(F.18), (F.22). The order £ and £ 2 terms are given 
by 

I'JT f "77" 

— A U{Qe + QruTB? + Z(Qe + Q m f B ) 2 } ~ ~. \Qe + T B Q m \ 2 & . (F.29) 

4T2B 1 ' 4t 2j b 

— * 

The linear term in £ do not give any contribution to the 1PI amplitudes. The contribution 
from the quadratic term gives 

\n(-^—\Q e + r B Q m \ 2 ) . (F.30) 

Thus the complete one loop contribution to the effective action is given by 

f !(r B ) = ln/( fc )(r B ) + In f (k \-f B ) + (k + 2) ln(2r 2B ) - ln(47r^ ) • (F.31) 

Up to an additive constant this agrees with the black hole entropy function for CHL 
models calculated in [27]. 

G The Integration Measure J(jb, £) 

— * 

In this appendix we shall compute the integration measure J{j B ,£) which arises from a 
change of variables from t±, r 2 to the normal coordinates: 

d 2 T 



(r 2 ) 2 



= J{f B ^)d%. (G.l) 



We first note that the duality invariant metric 

\(dr 2 + dr 2 ) (G.2) 

describes a metric of constant negative curvature —1. Since this is a homogeneous space, 
J7"(tb,£) cannot depend on t b . Furthermore with an appropriate change of variables we 
can bring the metric and measure to the standard form: 

d 2 r 

ds 1 = d9 2 + sinh 2 6d<p 2 , — = sinh0d0#. (G.3) 

T 2 
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Since J is independent of the base point tb, we can calculate it by taking the base point 
to be at 9 = 0. The geodesies passing through this point are constant (f> lines, and the 
length measured along such a geodesic is given by 9. Thus we have 



£ = (#cos0, 6* sine/)) . 



(G.4) 



This gives 



d 2 £ = 9d9d<P . 



(G.5) 



Comparing this with (G.3) we get 



d 2 r sinh 9 



d 2 i = i sinh ||| d 2 



(G.6) 



9 



Thus 



J(T B ,£) = ±rSmh\£\. 



(G.7) 
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